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Abstract: Supersymmetric extensions to the standard model provide viable dark matter 
candidates and can introduce additional charge-parity (CP) violation needed for obtaining 
the observed baryon asymmetry of the universe. We study the possibilities of scalar and 
neutralino dark matter with spontaneous CP violation in the next-to-minimal supersym- 
metric standard model with a right-handed neutrino. The observed relic density can be 
produced both by a neutralino or a right-handed sneutrino as the lightest supersymmetric 
particle but when CP is violated new annihilation channels become available and in general 
lower the relic density. We consider collider phenomenology for a number of benchmark 
points which all satisfy experimental constraints and have either the neutralino or the 
right-handed sneutrino contribute to the dark matter abundance. 
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1 Introduction 

One of the most compelling hints of physics beyond the standard model is the cosmological 
observation that more than 80 % of the mass of the universe is composed of dark matter 
(DM) [1]. This evidence comes both from the astronomical observations of gravitational 
effects at various scales and the cosmic microwave background measurements which are 
consistent with such a large amount of dark matter. 

Supersymmetry (SUSY) is one of the prime candidates of new physics and may provide 
an appropriate dark matter particle. Supersymmetric models generally contain terms, 
which violate baryon and lepton number conservation, and potentially lead to fast proton 
decay. The presence of such terms can be prevented by conservation of i?-parity [2-6]. 
The remarkable consequence of -R-parity conservation is that the lightest supersymmetric 
particle (LSP) will be absolutely stable. The LSPs left over after the Big Bang could then 
explain the observed dark matter relic density (RD). The LSP is thought to be uncharged 
(in both electric and color charges) so it interacts only weakly or through gravitational 
interactions. 

The candidates for the lightest supersymmetric particle in the minimal supersymmetric 
extension of the standard model (MSSM) spectrum are the lightest neutralino, gravitino, 
and sneutrino. The LEP-collider searches exclude a light left-handed sneutrinos as the LSP, 
and masses beyond LEP's reach are ruled out by direct detection searches [7-10]. However, 
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the superpartner of the right-handed neutrino is a viable dark matter candidate [11, 12]. 
A pure right-handed (RH) sneutrino has a very reduced coupling to ordinary matter due 
to the smallness of the neutrino Yukawa coupling and is therefore not produced thermally 
in the early universe. If, however, the right-handed sneutrino is by other means coupled to 
the rest of the observable sector it may be a thermal relic and provide for an appropriate 
dark matter relic density. 

Supersymmetry must be a broken symmetry at low energies. The explicit SUSY break- 
ing is introduced softly so that no quadratic divergences appear. This requires inclusion 
of all the possible gauge invariant breaking terms in the Lagrangian. The large number 
of breaking parameters parametrize the supersymmetry breaking, which is expected to be 
spontaneous in a more complete theory. Soft SUSY breaking introduces a large number 
of complex phases, which can lead to large CP violation effects [13]. While cosmological 
observations of the baryon asymmetry of the universe suggest that additional CP violation 
beyond what the Standard Model (SM) offers is necessary [14] , SUSY phases typically pro- 
duce excessively large electric dipole moments (EDMs) [15, 16]. One way to manage this 
problem is to impose CP conservation on the Lagrangian and introduce spontaneous CP 
violation (SCPV) [17]. SCPV is an attractive solution since it radically decreases the num- 
ber of free CP- violating parameters, in addition to allowing study of the interplay between 
different CP observables [18]. It is well known that SCPV is not possible in the MSSM 
and requires at least one and preferably two additional singlet fields [19, 20]. Therefore, 
we are led to consider SCPV in the next-to-minimal standard model (NMSSM). 

The NMSSM [21-23] provides also a solution to the so-called \i problem [24] by in- 
troducing an extra singlet scalar superfield S (lepton number L = 0), whose vacuum 
expectation value (vev) will generate an effective /i term. This naturally leads to fx to 
be of the order of the electroweak (EW) scale. The singlet also contributes to the Higgs 
masses already at the tree-level, which may lead to heavier Higgses than in the MSSM. 
Adding yet another singlet superfield N (L = 1) to the NMSSM provides for right-handed 
neutrino and sneutrino states. Since both S and N are gauge singlets it is tempting to 
have one field to do the job of both. This, however, leads to fine tuning problems as well 
as either explicit imparity violation or spontaneous lepton number violation with a super- 
fluous Goldstone boson in the spectrum [25]. Therefore, both singlets are assumed to be 
included in the model. Interestingly, a non-vanishing vev of the superfield S enables an 
effective Majorana mass term for the right-handed neutrino similar to the effective fi term. 
Moreover, the presence of the singlet S also leads to the electroweak scale interactions of 
the right-handed sneutrino with other MSSM fields. 

In this work we study neutralino and right-handed sneutrino dark matter in the frame- 
work of an extended NMSSM when CP is spontaneously violated. For the case of no CP 
violation, the possible right-handed sneutrino dark matter in the NMSSM has been earlier 
investigated in [26, 27] and neutralino dark matter in [28, 29]. Here we first introduce the 
changes in the model due to the spontaneous CP violation. We scan over interesting ranges 
for a number of parameters. After taking into account all relevant experimental constraints 
it is seen that CP violation decreases the dark matter relic density because of new available 
decay channels. The possible collider signals are considered for several benchmark points, 
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two of which have a sneutrino as an LSP and five of which have a neutralino LSP. 



2 NMSSM with a singlet neutrino superfield 

We start by introducing the important properties of the model. The superpotential for the 
NMSSM with singlet superfields S (L=0) and N (L=l) is given by 

W = e a/3 (YgHZtfEj + Y^H^D 3 + YffQfH^ + Y^LfH^N + \SH%H? 

+X N NNS+^SSS, (2.1) 

where e a p (a, /3 = 1, 2) is a totally antisymmetric tensor with e\2 = 1, X denotes a superfield 
and the rest are dimensionless couplings. The soft SUSY breaking terms are 

V soit = [e af} (A%Y$WffiE s + A%YZH?Qlb. + A^YffQfH^ + A^LfH^N 
+ A X XSH^H^ + A Xn X n NNS + ^SSS + fs} + h.c. 
+ .Ui,,<l>;<I> ; + Ml^-e.e* + m\ x H\H x + m\ 2 H\H 2 + m 2 s SS*, (2.2) 

where 3> = {L, Q},@ = {E,N,U,D} are the scalar components of the corresponding 
superfields. A discussion on the theoretical merits of this and similar models (specifically 
the soft S tadpole) can be found in [23, 30, 31]. For phenomenological merits of this and 
similar models (NMSSM + right-handed neutrino) see [32-34]. 

In this work we take only the third generation Yukawa couplings to be non-zero for 
Yu/d/e an d impose at the electroweak scale Y^ = Y^ and A % N = An for all three genera- 
tions. We also take all soft scalar masses diagonal, Mf- = M 2 <5y, and take Mjj = = 
Mq, as well as Me = Ml- It is worth pointing out a tension in choosing natural sizes 
for some of theses parameters. If we want the \i problem solved then (S) ~ EW scale 
and it follows that M^ ~ (S) ~ EW scale, assuming that A and Atv are 0(1). Now the 
left-handed neutrino masses are generated by the conventional seesaw mechanism: 

M » = ( °t Tf ) m "L - -m D M^ l m T D , (2.3) 
y m x D M N J 

where mp ~ Yj^x EW scale and thus m UL ~ Y^x EW scale, meaning we need to choose 
Y N ~ 10~ 6 . Otherwise, even for (S) ~ Mjy ~ Mpi anc ^ we would have Yn ~ 10 2 on top 
of which we would have to choose A ~ 10~ 17 . 

Spontaneous CP violation is introduced by complex vacuum expectation values of the 
S and H® fields, 



W=L^J, (S)=v S e^. (2.4) 





A complex phase of the vev of H® can be absorbed by field redefinitions ({H®} = v%). 
Non-zero vevs for left- or right-handed sneutrinos would introduce spontaneous R-parity 
breaking. Since we are interested in viable dark matter candidates we consider only PL- 
parity conserving models. Thus there are only two complex phases 5s and 82 which enter 
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couplings and mass matrices wherever the vevs of S and appear. In particular the mass 
matrices of neutral scalars are no longer block diagonal with respect to the division of CP 
even and CP odd gauge eigenstates. The sneutrinos form an 8 x 8 mass matrix divided 
into 4x4 submatrices 

/ 2 2 \ 

rut* m 



M i = 2 T (2-5) 



where the subscript denotes CP-even/odd states. We introduce the following notation 





= Y^(A l N V2 cos 82 ± 2\nV2Vs cos(<52 — 8s) ± \v1v2 cos 8s), 


(2.6) 




= Af^ 1 (cos ->■ sin) , 


(2.7) 


c ±± 


= 2A^^AA iV ^'s , cos 8s ± 2kXnv 2 cos 25s ± 2AAati>i?;2 cos fe, 


(2.8) 




= C ±=t (cos ->■ sin), 


(2.9) 


m L,ij 


= Ml tij + Y£Y£(w? + + \m\ cos 2/3%, 


(2.10) 


m R 


= m^ + ^yaV 2 + ^) + 4A^I. 

i 


(2.11) 




The 4x4 submatrices are then 

ml = ( ^ ^ ) , mL = c+ _ ) , (2.12) 

"4 = I ^ P + j, rnl = (mlf. (2.13) 

where i,j are the family indices of the left-handed sneutrinos. The mixing between left- 
and right-handed sneutrino states is suppressed by the smallness of neutrino Yukawa Yv 
in and . The CP- violating mixing between even and odd states by can be 
further suppressed if the CP-violating phases are small. There is no mixing between even 
and odd left-handed states and the mass splitting for left-handed states is proportional to 
Yjy. Consequently our model has almost purely left-handed mass degenerate sneutrinos 
which are CP eigenstates to a high degree of accuracy. 

The dominantly right-handed sneutrinos are split by and D h into 

mj 12 ~m 2 R T VC+- 2 + D~+ 2 . (2.14) 

From Eqs. (2.8) and (2.9) we see that the mass splitting is proportional to A at but unsup- 
pressed by possibly small CP phases due o the complimentary nature of C" 1 and D 
with respect to the CP-violating phases. The mixing between even and odd states is done 
by D which may be suppressed by small CP violating phases. 

The Higgs mass matrix is also no longer block diagonal but becomes a 6 x 6 matrix 
where the CP-even and CP-odd states mix. The usual vacuum stability conditions are 
solved for the EW scale mass parameters m^, m// 2 , and ms- The addition of CP phases 
introduces new vacuum stability conditions which we solve by fixing the soft trilinear terms 
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A\ and A K . At tree level 

dV sin(^ - 25s) fl) 

= ^ A x = -Kv s . . . , (2.15 

ad 2 sin(d 2 + <5s) 

<9F 3KAuiV2V5sin(5 2 - 25 s ) + £ 3 sm5 s , 

— =0^A K = 2 . q c • (2.16) 

In our numerical calculations we use the 1-loop effective scalar potential including the 
corrections from the third generation (s)quarks to derive both the mass matrices and the 
vacuum stability conditions. 

The neutralino mass matrix is of the conventional NMSSM form with the addition of 
complex phases 

/Mi -am ffi^ e -«a \ 



M x o 



Mi 








M 2 


givi 


92«1 




V2 

















9i»l 
92 fl 




91^2 „-i<5 2 
92«2 „—■ i5 2 

-Xv s e l5s 





-At^ 2 







— Avi 




— Aw i 


2ra s e i<5s 



(2.17) 



V 

3 Constraints and the parameter space 

3.1 Parameters 

We carry out our numerical analysis of the effect of the CP-violating phases by randomly 
sampling the parameters that affect the Higgs and sneutrino spectra as well as the interac- 
tions of the dark matter candidate. The gaugino masses were left fixed to better illustrate 
the effect of bino dominance in the neutralino LSP. The parameters affecting squark masses 
were chosen such that we avoid the current experimental limits completely. We have gen- 
erated two data sets, one CP-conserving {62 = 5s = 0) and one where we vary the CP 
phases along with the other sampled parameters. Table 1 gives the values for the fixed 
parameters and the ranges for the sampled parameters. 

3.2 Tools 

The computational tools we use are LanHEP [35] for model file creation and micrOmegas 
(v. 2. 4.1) for the main numerical analysis [36]. Calculating the spectrum and diagonaliza- 
tions are performed using the EISPACK [37] routines and calculation of electric dipole 
moments and rare decay branching ratios is performed by a dedicated Mathematica code. 
B-physics constraints are calculated by NMSSMtools [38, 39] using the supplied NMSSM 
model. 

There are two caveats to the use of NMSSMtools which we will briefly address here. 
One is that the NMSSM model does not include a right-handed neutrino. However, loop 
contributions to the calculation of BR(i? — > X s ~f) and BR(i? s — > /i + fi~) that involve 
a right-handed sneutrino will be suppressed by Y^/Y^, compared to the corresponding 
left-handed contributions and are thus negligible in our model. The other caveat to using 
NMSSMtools is that CP violation is not implemented in the supplied NMSSM model. Two 
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Mi 


300 GeV 


M 2 


600 GeV 


M 3 


1800 GeV 


Mq 


1000 GeV 


A t 


1500 GeV 


A b 


1500 GeV 


A T 


-2500 GeV 




10~ 6 


A Nl 


GeV 



tan /3 


2-50 


1* 


- 500 GeV 


X 


- 0.8 


K 


- 0.8 


A 

A\ 


-i i~\t~\r\ /~1 T T -t r\r\r\ /"I t t 

-1000 GeV - 1000 GeV 


A 

A R 


-1UUU Ljev - lUUU uev 


v s 




An 


- 0.8 


A 


-100U Lrev - 100U Lrev 


M N 


- 500 GeV 


M L , E 


- 500 GeV 


8 S 


- 2vr* 


52 


- 2vr* 




-1000 GeV - 1000 GeV 



Table 1. The parameters and sampling ranges used in our two data sets. *EDM considerations 
limit the range of the phases considerably and we use the ranges [0,0.3], [ir — 0.3, n + 0.1], and 
[2n — 0.1, 2n] for generating points that are subject to all of the experimental constraints. 



factors can affect the calculation: changes in the model's spectrum due to non-zero CP 
phases, and changes in the couplings of the fields. For the particles appearing in the tree 
level and 1-loop contributions to BR(B — > X s j), BR(B + — > t + u t ), and BR(£? S — > 
the spectrum changes continuously from #2, Sg = — )■ 62, 63 ^ 0. To judge the impact of CP 
phases on the calculations we did a comparative analysis using the MSSM for which both 
CP-conserving and CP-violating code bases exist in NMSSMtools. In the limit of small 
CP phases we find that the difference in the results for BR(i3 — > X s j), BR(i3 + — > t + v t ), 
and BR(5 S — > is at the 5% level for phases as large as 7r/10 and more generally at 

the level of 1% for phases < 0.1. 

This gives us confidence that the CP-conserving NMSSM model can be used to cal- 
culate -B-physics constraints in the regime of small CP-violating phases (82,63 < 0.1). 
Satisfying the EDM constraints restricts us naturally to this regime. We find that in our 
CP-violating data set most points that satisfy EDM constraints have l^l < 0.05 (modulo 
7r). The phase 5s is not as constrained and in general the phases 82 and 83 can both have 
large values such that the physical phases entering EDM calculations are fine tuned to 
be small. However, to simplify calculations and to be able to use the NMSSMtools pack- 
age mentioned above we restrict both phases to the ranges [0,0.3], [ir — 0.3, tt + 0.1], and 
[2tt — 0.1, 2ir] when generating points that are subjected to all the experimental constraints 
we impose. 

3.3 Constraints 

In order to examine only physically relevant parameter regions, a number of constraints 
are taken into account during our computations. We check the scalar sector for vacuum 
stability and the sparticle spectrum against the most recent PDG limits [40]. 
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The LEP experiment places the lower limit on the standard model Higgs boson mass 
at 114 GeV. As the Higgs boson can only be detected by its decay products, not directly, 
a model independent way to use the experimental results is to investigate the bounds on 
Higgs boson couplings to ZZ and hZ. In the NMSSM the couplings hZZ and hhZ may 
differ from the standard model and consequently the mass of the Higgs boson can be below 
114 GeV, provided that the couplings are sufficiently small. The reduced Z couplings are 

9 hi zz = cos p On + sin p 2i (3.1) 
g hihjZ = (sin P 4i + cos P 5i )(cos p0 2j - sin pO^-ii^j) (3.2) 

with Oij the 6x6 neutral Higgs mixing matrix. We use the most recent limits these 
couplings impose on the Higgs mass, which are presented in [41]. The recent results from 
the LHC and Tevatron experiments imply that the SM-like Higgs boson mass should be 
around mu ~ 125 GeV [42, 43]. In our analysis we thus also require one of the Higgs 
states to have a mass between 123 - 128 GeV. The 125 GeV Higgs signal in the NMSSM 
has recently been studied in Refs. [44-47]. 

The phases 5 2 and 5s that introduce CP violation in our model are flavor diagonal and 
thus are expected to be highly constrained from, e.g, the electron electric dipole moment 
measurements. Current limits [48] constrain the electron EDM to be below 

d e < 1.05 • 10~ 27 ecm. (3.3) 

The NMSSM specific contributions to BR(i? — > X^) from the extended Higgs and neu- 
tralino sectors arise only at two-loop level [23] and, in general, are considered to be small 
[49, 50]. In our calculations, the theoretical uncertainties [51, 52] are combined with the 
present experimental value [53], which leads to [54] 

BR(B ^ X sl ) = (355 ±142) x 10~ 6 . (3.4) 

The impact of the NMSSM on the branching ratio BR(Z? + — > t + v t ) is only indirect, 
and the extra singlet contributions are suppressed by the small neutrino Yukawas as men- 
tioned earlier. The new physics contribution to the branching ratio can be quantified by 
defining a ratio [55, 56] 

n on ^ BR(B+ -» T+t/ r )sM+SUSY q -| q ,„ -s 

0.99 < pp , p+ - , — v < 3.19, (3.5) 

where the numerator denotes the branching ratio in the SUSY scenario. The constraint 
(3.5) tends to prefer small values of tan P in order not to decrease the ratio too much below 
the lower limit. A large charged Higgs mass decreases the new physics contributions in 
general. 

The effect of our model on BR(S S —> /i + /i~) at one loop comes from right-handed 
(s) neutrinos and is thus suppressed by the smallness of the neutrino Yukawa coupling. We 
apply the limit from [57] 

BR(B S -> < 4.5 • 10~ 9 (3.6) 
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The anomalous magnetic moment of the muon has been measured quite precisely. How- 
ever, there is uncertainty in the reliability of the theoretical prediction due to hadronic and 
non-perturbative effects. Therefore, we do not use the magnetic moment as a constraint. 

The relic density of cold dark matter in the universe is determined by the Wilkinson 
Microwave Anisotropy Probe (WMAP) [1] to be Q c h 2 = 0.1126±0.0036. If 10 % theoretical 
uncertainty is added [58] we find the preferred WMAP range 

0.0941 < n c h 2 < 0.1311 (3.7) 

at 2 a level. As the dark matter may also contain a component of a non-supersymmetric 
origin, we have used only the upper bound as a constraint. 

4 The effects of SCPV on the relic density and other constraints 

In this section we show the effects of CP violation on the relic density of the LSP. All plots 
contain points from the CP-violating data set unless otherwise specified. All points also 
satisfy PDG constraints on the mass spectrum and vacuum stability. 

Spontaneous CP violation is known to change the spectrum of the scalar sector in a 
discontinuous manner. For small values of the phases there appears a light state, hs, in 
the spectrum [59] . Experimental constraints on the Higgs boson mass force this light state 
to be singlet dominated [41]. We define the doublet content of hs as 

e sd = ^ Of 2 (4.1) 
i=l,2,4,5 

with O the neutral Higgs mixing matrix and the free index i runs over the elements cor- 
responding to the lightest Higgs' gauge eigenstates separated into CP-even and CP-odd 
parts, i.e. {Reii^, ReH® , ReS, ImH® , ImH® , ImS}. We find that for our data points the 
value of e sc i is between 0.1 - 10~ 5 . The presence of this light singlet state has different 
effects depending on what kind of LSP we have. 

The addition of the singlet to the model also opens up the possibility for the neutralino 
LSP to have a significant singlino admixture and thus behave differently from its usual 
MSSM counterpart. We denote the singlino component of the neutralino LSP as 

e xs = l^sil 2 , (4.2) 

where ./V is the neutralino mixing matrix. We find that already at values of 0.1 the be- 
havior of, e.g., the relic density of the neutralino depends significantly on the Lagrangian 
parameters associated with the singlino. 

The other masses in the spectrum are also affected by the CP phases, but in a con- 
tinuous way, meaning that for small values of the phases the spectrum is very close to the 
CP conserving spectrum. This means that in calculations of observables, such as the relic 
density or rare decay branching ratios, the dominant effect comes from the new light scalar 
introduced by SCPV in our model. 
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We already discussed that for B physics constraints the effect of small CP-violating 
phases is negligible. We additionally find that the appearance of hs has very little ef- 
fect since it doesn't appear in the calculations of BR(£? — > Xg'y), BR(i3 + — > t + u t ), and 
BR(S S — > until two loop order. The charged Higgses contribute already at lower 

order but their masses change continuously with the CP phases and thus the effect of CP 
violation remains limited in the regime of small CP phases. However, in the case of the 
relic density there are new, possibly dominant, annihilation channels that open up with 
the appearance of hs depending on what kind of LSP we have. 

Figures 1 and 3 show the leading coefficients for two classes of direct annihilation 
channels, i.e., those with //, / € [q, I, u], and those with hshs final states in the case of a 
neutralino or a sneutrino LSP, respectively. We show neither the diagrams for vector boson 
final states nor for the various co-annihilation channels; the coefficients are all sub leading 
or the same as for the // diagrams. In the numerical calculations we naturally take all 
these channels into account. As we can show, this rough approximation provides a good 
understanding of the parametric dependence of the relic density. We have abbreviated the 
typical electroweak coupling ttt-^ — = C ~ 0.37. 




~c 2 



~ 



d 



d 



C-(Y f ,X N e sd ) 
k • (\N,Y f e sd ) 




~C 2 e 2 



sd 



~ Ce sd K 

,2 



d 



d 



(\,Ce sd ) 2 
(K,Xe sd ) 2 



Figure 1. Tree level pair annihilation channels for the neutralino into // /hshs- Leading coeffi- 
cients for the doublet (d) and singlet (s) component of the neutralino are given. The addition of a 
light singlet scalar to the spectrum opens additional channels. 



Fig. 1 shows the pair annihilation channels in the case of a neutralino LSP. We have 
separated the leading coefficients for the case of a pure singlino as "singlet" and group the 
couplings to all the other gauge eigenstates under "doublet" . We see that for the doublet 
component of the neutralino the dominant contributions are ~ C 2 with most of the new 
hs final states suppressed by e sd . 

For the singlet component of the neutralino the story is different. The unsuppressed 
~ Aat, channels for the // final states are available only if the right-handed neutrino 
is lighter than the LSP. Even then, m UR depends on Xjy directly so having m UR small may 
also significantly suppress these channels if, e.g., {S) is large. The channels for hs final 
states on the other hand contain unsuppressed k 2 dependence. However, k also controls 
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the neutralino mixing in a way that for small values of k the LSP becomes highly singlino 
dominated. In that case, the only available annihilation channels are those now suppressed 
by the smallness of k, which leads to a larger relic density. 

These effects are shown in Fig. 2. In the left plot we can see that for small values 
of k < 0.05 none of the neutralinos have an appreciable doublet component and also that 
the relic density for singlino dominated LSPs is controlled by k. In the right plot we see 
how with a growing singlino component the presence of the new hs final state channels 
significantly decrease the relic density compared to the CP-conserving scenario beginning 
a t e x s ~ 0.001. We also see that once only the singlino channels are available the relic 
density can again become very large due to suppression by k. 




0.01 0.02 0.05 0.10 0.20 0.50 1.00 10" 4 0.001 0.01 0.1 1 



Figure 2. Left: The relic density against the trilinear coupling k for neutralino LSPs. Right: 
The relic density against the singlino component of the neutralino LSP. Points with (blue circles) 
and without (orange triangles) significant singlino component. Green boxes depict points from the 
CP-conserving data set. The grey band indicates the the current WMAP limits on the relic density. 

Fig. 3 shows the pair annihilation channels for the sneutrino LSP. In the case of a 
left-handed sneutrino the channels decaying into hshs are suppressed by additional factors 
of C 2 and/or e s d- It is known that the left-handed sneutrino is a poor LSP candidate if one 
wants to saturate the dark matter relic density limit [60]. This is caused by the fact that 
the dominant ~ C 2 weak decay channel is always present for the left-handed sneutrino. 

In the case of a right-handed sneutrino Fig. 3 shows that the hshs channels are all 
of the same order as the // channels and collectively depend to leading order on Atv or 
A 2 ^. Thus, it is expected that the relic density will be lower when compared to the CP- 
conserving case. The left plot in Fig. 4 shows how the relic density for right-handed 
sneutrinos depends strongly on An, while the relic density for left-handed sneutrinos stays 
almost a constant. The right plot in Fig. 4 shows the drop of the relic density for the 
CP-violating case for right-handed sneutrino LSPs when compared to the CP-conserving 
data. Although much of the set overlaps there is a clear systematic drop in the mean relic 
density for the CP-violating set. 

In Fig. 5 all of the aforementioned effects are shown in unison. In the plot on the 
left the fixed gaugino mass M\ = 300 GeV causes the neutralino relic density to spike at 
m LSP = M± due to the small bino annihilation process cross section. We also see that for 
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v y x h s 



Figure 3. Tree level pair annihilation channels for the sneutrino into ff/hghs- Leading coefficients 
for the left-handed (L) and right-handed (R) case are given. In the left-handed case the addition of a 
light singlet to the spectrum opens up new channels but they are all suppressed by additional factors 
of C 2 or e 2 ad when compared to the originally dominant weak decay channel. In the right-handed 
case the new hs final state channels are of the same order as the other decay channels. 



a 




Figure 4. Left: The relic density against the trilinear coupling Ajv for right-handed (blue circles) 
and left-handed (orange triangles) sneutrino LSPs. The grey band indicates the current WMAP 
limits on the relic density. Right: A box-and- whisker diagram of the relic density against the trilin- 
ear coupling A n comparing the distribution of right-handed sneutrino LSPs in the CP-conserving 
(light green) and CP-violating (blue) cases. Boxes are the 25 - 75 percentile range and whiskers 
denote the complete range of the data. 



other values of m LSP the neutralino has to be singlino dominated (yellow dot) if we want 
to saturate the relic density limit and how in the case of CP violation the relic density 
is significantly reduced compared to the CP conserving case. In the right plot we see 
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Figure 5. The relic density against the LSP mass. Left: for neutralino LSPs in the CP-conserving 
(light green triangle), CP- violating doublet dominated (orange triangle) and CP- violating significant 
singlino admixture, i.e. e xs > 0.1, (yellow dot) cases. Right: for sneutrino LSPs in the CP-violating 
left-handed (orange dot) and CP-violating right-handed (yellow triangle) case. The CP-conserving 
left-handed (light green dot) set overlaps completely with the CP-violating counterpart. In both 
plots all points satisfy PDG constraints on the mass spectrum and vacuum stability, large brown 
triangles or circles indicate points that pass all of the other experimental constraints we impose. 
The grey band indicates the current WMAP limits on the relic density. 



that left-handed sneutrino LSPs are poor candidates for saturating the relic density limit 
regardless of CP violation, whereas for right-handed sneutrinos this is not a problem. The 
scattered very low relic density points for the left-handed sneutrinos are due to the presence 
of efficient slepton co-annihilations. Both plots also show how few points pass all of the 
experimental constraints we impose. We find that the most restrictive constraint is the 
electron EDM, e.g., in Fig. 5 removing all other CP, B physics, and rare decay constraints 
only yields a single additional "allowed" point. 

In general we find that for constrained phases the effect on the relic density from 
variations in the spectrum of the model is much smaller than the effect introduced by 
the appearance of a light singlet scalar providing new annihilation channels for the LSP. 
Depending on the LSP and its composition these channels are governed by the Lagrangian 
parameters n and Xn- 



5 Direct dark matter searches 

There are currently several ongoing searches of weakly interacting massive particles ( WIMPs) 
in direct detection experiments. The goal of the experiments is to measure the WIMP 
nucleus interaction in the detector material. WIMPs scatter when they interact with nuclei, 
and the recoil energy can then be measured. As yet, the XENON100 experiment puts the 
tightest limits on the WIMP-nucleon scattering cross section [61]. To compare the cal- 
culated proton/neutron cross sections with the experimental limits, we use a normalized 
cross section for a point-like nucleus [62] : 

(z^r p + (A-z)^) 2 

"SI = ^ ' t 5 - 1 ) 
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Z and A being the atomic and mass number of the target element and <r njP the spin- 
independent cross-sections for neutron and proton target, respectively. Because there are 
large uncertainties in the local density of dark matter and in the nuclear matrix elements 
that enter the computation of <tsi [63, 64], the direct detection limits are only indicative. 

Fig. 6 shows the LSP mass versus the spin independent cross section for a set of CP- 
violating points, as well as the XENON100 limit. We see that the XENON100 limit actually 
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Figure 6. Spin-independent WIMP-nucleon cross section for CPV points. Points with triangular 
shape (black) have the right-handed sneutrino as an LSP, while the dots (orange) have neutralino 
LSP. The (green) cross above a point corresponds to a relic density above the WMAP upper limit, 
hence an excluded point. The solid line represents the XENON100 experimental limit. 



serves as a constraint, since a relatively large portion of the parameter points are ruled 
out by the direct detection search. A roughly equal number of points are still allowed by 
the XENON100 limit. The WMAP upper limit constraint removes most of the sneutrino 
LSP points allowed by the XENON100 limit. However, some points still evade both dark 
matter constraints. A large number of neutralino LSP points are still allowed in the low 
neutralino mass region. In ref. [65] it was found that in the CP-conserving NMSSM the 
direct detection limits for the neutralino LSP still allow a large number of points with low 
fine tuning, especially for a large A parameter. 



6 Signatures at the colliders 

In this section we focus on the possible signals from the supersymmetric particle spectrum 
realized by the model respecting all relevant low energy constraints and generating an 
acceptable amount of relic density. Notwithstanding the fact that we have already restricted 
the parameter space of the model significantly by various constraints, recent observation 
of a ~ 125 - 126 GeV Higgs like boson at LHC by both the CMS [66] and ATLAS [67] 
Collaborations further helps in narrowing down the allowed sparticle spectrum. 

^or Xenon: A = 131, Z = 54. 
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As discussed in Sec. 4 in the model used here with SCPV a light scalar state hi = hs 
exists and thus the branching ratio of the decay of the ~ 125 GeV Higgs, h>2, into a pair hs 
is of particular interest. A limit on BR(/i2 — > hshs) could be useful in further constraining 
this type model. We find that a priori requiring mh s > is problematic since a larger 

mass for the singlet type Higgs implies larger CP violation and thus the limits on the EDMs 
exclude these points. Instead we extract the hhh coupling as 



with O the Higgs mixing matrix and require \C2x1 \ < lGeV. From the points that clear 
this limit we have selected a few which also clear the EDM, B-physics, and rare decay 
constraints and are favored by relic density and analyze their signal at LHC. We have listed 
the benchmark points (BP) in Table 2. Note that LHC has already put some strong limits 
on the squark and gluino masses [68, 69], albeit these constraints are model dependent 
and also depend on the way the heavy particles decay. However, for our analysis, we 
have made a consistent choice of parameters where the strongly interacting sector of the 
sparticle spectrum remains above a TeV and is, therefore, unconstrained from LHC data. 
This implies that our analysis will focus on the weakly interacting sector of the spectrum 
viz., sleptons, gauginos/higgsinos and the Higgs. 

Our choice for parameter points leads to a right-handed sneutrino LSP or a neutralino 
LSP. For BP6 and BP7 the LSP comes out to be a right-handed sneutrino with a mass of 
~ 6 GeV and ~ 184 GeV, respectively. All the other benchmark points have a neutralino 
as an LSP. We focus on each benchmark point individually and highlight the most likely 
signals for the different parameter choices at LHC in our model. 




(6.1) 



<j) e {ReH",Rei^,ReS,ImiTi,Im#£ 



ImS} 
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BPl 


BP2 


BP3 


BP4 


BP5 


BP6 


BP7 


tan p 


45.7 


12.4 


33.9 


26.9 


30.1 


43.5 


44.9 


A 


0.114 


0.355 


0.363 


0.347 


0.430 


0.259 


0.341 




0.038 


0.63 


0.231 


0.294 


0.357 


0.174 


0.292 


vq (GeV) 


4277.2 


1160.2 


964.8 


914.8 


487.5 


1881.9 


816.7 




0.721 


0.008 


0.032 


0.498 


0.269 


0.035 


0.750 


A\ M (GeV) 


337.0 


25.2 


-668.9 


319.7 


365.3 


-135.5 


-975.5 


Mat (GeV) 


447.7 


449.8 


494.7 


401.1 


341.5 


207.6 


135.7 


M l ,b (GeV) 


307.1 


419.5 


432.0 


431.9 


472.7 


320.3 


377.0 


*s 


3.228 


0.156 


3.142 


3.203 


3.213 


3.199 


3.173 




0.111 


0.034 


0.010 


0.249 


3.037 


0.142 


0.173 


£ (GeV) 


118.1 


-991.3 


522.9 


318.9 


253.0 


109.8 


300.6 


LSP-type 


Xo 


Xo 


Xo 


Xo 


Xo 






% SP (GeV) 


296.3 


290.9 


286.3 


276.7 


194.6 


5.7 


183.7 


n c h 2 


0.062 


0.098 


8.3- 10" 3 


8.3 ■ 10~ 4 


6.5 • 10~ 4 


3.3 • 10" 6 


2.9 ■ 10" 4 



Table 2. Benchmark points for studying the signals of the supersymmetric particles at LHC. 
Note that fi = \v$ and other parameters are the same as shown in Table 1. 



• BPl 

This point gives a neutraiino LSP with mass ~ 296 GeV while the lightest sneutrino is 
nearly degenerate at ~ 300 GeV. The chargino xf 1S much heavier (470 GeV) than the 
second neutraiino \\ (322 GeV). Although this point is allowed by all experimental 
data it does not give any significant cross sections at LHC in any channel. The 
chargino pair production which is the largest is about 8 fb at LHC with yfs = 14 
TeV. The dominant decay for the chargino is to charged leptons and lighter sneutrinos 
which then decay invisibly. So the final state which could be observed at LHC is 
£f£j$T which would require substantial luminosity to separate it from the large SM 
background coming from the W + W~ production. The point also gives a very light 
scalar of ~ 2.7 GeV which is dominantly singlet with very suppressed couplings to 
all SM particles, such that it is not constrained by any experimental data. 

• BP2 

For this point we have a relatively small tan f3 and Ajv as compared to all other points. 
The LSP is the lightest neutraiino which has a mass of ~ 291 GeV. The small value 
of \n leads to a very light right-handed neutrino of mass ~ 19 GeV. This state has a 
very small mixing with the left-handed neutrino and thus is not constrained by the 
measured Z width or other low energy data. The strongest signal for this point at 
LHC is through the production of first chargino and second neutraiino {X1X2) w h° se 
mass is about 400 GeV with a cross section of ~ 5 fb for the current run at LHC 
with yfs = 8 TeV while the next significant production cross section is for the pair 
production of charginos (xiXi) with a cross section of ~ 4 fb. All other production 
modes are suppressed. The cross sections are increased by a factor of ~ 3 — 4 at LHC 
with yfs = 14 TeV while other modes also become accessible with high luminosity. 
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However, we note that the mass difference of these states with the LSP is only about 
100 GeV. Both the xt an d X2 decay to the LSP and a weak gauge boson with 
100% branching probability. Thus the pair production of charginos will lead to a 
final state of W + W~ $ T where both W's will have small transverse momenta. The 
X1X2 production leads to final states with W + Z$ T with no suppression in branching 
ratios. This leads to interesting signal of a Z peak, an associated charged lepton and 
large missing energy. The SM background will be mostly from WZ and triple gauge 
boson production where the large WZ production can be suppressed by demanding 
large 

• BP3 

This point gives a neutralino LSP with a mass of ~ 286 GeV with next lightest 
supersymmetric particles being \\ (338 GeV) and xt (341 GeV). This point also 
gives us a light right-handed neutrino with mass 61 GeV and an ultralight scalar 
with mass 300 MeV. The largest production cross section in this case is again xtx\ 
and X1X2 but the xt decays via off-shell weak gauge boson (W + ) and off-shell 
sleptons/sneutrinos leading to final state particles with small transverse momenta. 
The x% decays to an ultralight singlet scalar and the LSP with 100% probability. The 
signal production cross sections are about ~ 7.5 fb and ~ 3 fb, respectively, at LHC 
with y/s = 8 TeV, and a factor ~ 3 - 4 larger at LHC with y/s = 14 TeV. As the final 
decay products do not carry large transverse momenta, isolating the signal from the 
background will be very challenging. 

• BP4 

This point has a very similar spectrum to BP3 except for a much heavier right-handed 
neutrino with a mass of ~ 910 GeV as A at is much larger. The lightest scalar is the 
dominantly singlet component of the Higgs with mass ~ 17 GeV. The neutralino LSP 
has a mass of ~ 277 GeV. The next to lightest supersymmetric particle (NLSP) is 
the xt (3H GeV) while x% (316 GeV) also is close and nearly degenerate. The mass 
difference again forces the decay similar to that for BP3 with final states with small 
transverse momenta. The largest cross section in this case is for the pair production 
of charginos (xtxi) which is about 13 fb at LHC with y/s = 8 TeV which increases 
to ~ 40 fb for y/s = U TeV. 

• BP5 

This point corresponds to a much lighter supersymmetric spectrum with a neutralino 
LSP of mass ~ 195 GeV. The NLSP is X \ (1" GeV ) with Xi ( 204 - 5 GeV ) the next 
lightest. It also has a 262 GeV right-handed neutrino. The light spectrum leads to 
much larger production cross sections, and the dominant production channels in this 
case are XiXi^xtxi ail d xtx2 with production cross sections of around 74 fb, 86 
fb and 82 fb, respectively, at LHC with y/s = 8 TeV. However, the mass difference 
of xt an d X2 wr th Xi i s l ess than 10 GeV, and thus the transverse energy of the 
decay products is small. It would be quite difficult to trigger on such soft leptons 
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and jets that result from the 3-body decays of X\ an d xi- A detailed analysis would 
be required to isolate signals which would be relevant to highlight for LHC searches. 

• BP6 

Unlike the other points, we get a right sneutrino LSP for this parameter choice with a 
mass of ~ 6 GeV. Pair production of the LSP has the largest cross section but will pass 
through the detector undetected. The other significant cross sections (xi~Xi~i X1X2) 
are relevant only at LHC with ^/s = 14 TeV with production cross sections of less 
than 8 fb. The decay of the Xi aR d X% is however different from that of BP2 with 
suppressed probabilities in similar channels and therefore the signal gets suppressed 
further. Much like BP1 we get a light scalar with a mass of ~ 2.5 GeV. We also get 
a light right-handed neutrino of mass 130 GeV which can lead to interesting signals 
if there is significant mixing with the left-handed neutrinos. 

• BP7 

Our final benchmark point gives us a right sneutrino LSP of mass ~ 184 GeV. This 
point also gives a light scalar Higgs of mass ~ 8.4 GeV in the mass spectrum. The 
Xi (254 GeV) is the NLSP while the right-handed neutrino mass is above 1 TeV. The 
largest production cross section in this case is for X\X\ where the chargino mass 
is ~ 273 GeV, which at LHC with ^/s = 8 TeV is about 22 fb and jumps up to 
~ 65 fb for -y/s = 14 TeV. The chargino decays to a r and LSP with 60% branching 
probability while it decays 20% of the time to electron/muon and LSP. This leads to 
final states of two charged leptons with large missing energy. Other relevant channels 
are from production of X1X2' Xi Xi an d X^X2 where the Xi decays invisibly while the 
X2 decays to the light scalar Higgs and Xi with 100% probability. 

Although we discuss some possible signals at LHC for our model with different choices 
of the parameter space a more detailed analysis of both signal and background is needed 
which we leave for future work. The CP- violating effects due to non-zero phases can affect 
several observables reconstructed from the cascade decays of the charginos, neutralinos or 
the scalars in the model and have been highlighted in various studies [70-73] and we aim 
to apply these studies in our future analysis of the model at LHC in our next publication 
[74]. 

7 Summary and discussion 

We have studied the viability of the right-handed sneutrino and neutralino dark matter 
in the NMSSM, and the effect of CP violation on the dark matter relic density and the 
direct detection potential. The studied model contains, in addition to the MSSM fields, 
two singlet superfields and interaction terms. This extends the Higgs sector, adds the 
right-handed neutrino field to the model, and allows the right-handed sneutrino to become 
a thermally produced dark matter candidate. We took into account the BR(1? — > Xsj), 
BR(Z? + —J- t + v t ), and BR(£> S — > bounds, the discovery of the Higgs and sparticle 
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searches, and constrained the model parameters by calculating the dark matter relic density, 
including co-annihilations. We calculated also the direct detection rates for the nucleon- 
DM spin-independent cross section, and compared that to the XENON100 limit. We have 
found that for constrained CP- violating phases the effect on the relic density from variations 
in the spectrum of the model is much smaller than the effect introduced by the appearance 
of a light singlet scalar providing new annihilation channels for the LSP. In particular for 
neutralinos and right-handed sneutrinos these new channels lower the relic density. In some 
cases, the left-handed sneutrinos co-annihilate very efficiently with other sleptons resulting 
in very low relic density. The spin-independent cross section can be used to constrain the 
models with a light sneutrino LSP, as for these models the direct detection limit tends to 
be complementary to the relic density constraint. The neutralino dark matter, however, 
remains rather unconstrained by the direct detection limits. We have also used our analysis 
to identify a set of benchmark points which satisfy all constraints and favor DM data. We 
have used these benchmark points to highlight the most important signals for the model 
that can be observed at the current and future run of the LHC. 
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